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SECTION  I 


INTRODUCTION 


An  important  factor  in  the  execution  of  military  ope.rations  is  the 
weather  at  the  scene  of  the  battle.  Whether  or  not  air  operations  can 
be  undertaken  (and  which  ones)  depends  on  the  ceiling  and  visibility 
conditions.  Therefore,  many  weapons  allocation  models  have  a  capability 
for  simulating  these  conditions.  These  conditions  are  determined  for 
the  most  part  by  launch  envelopes  of  aircraft-missile  systems.  For 
Instance,  steep  dives  require  high  altitudes.  If  the  ceiling  or  visibility 
is  low  so  that  the  pilot  cannot  see  the  target,  then  the  maneuver  cannot 
be  accomplished.  It  is  therefore  desirable  to  have  a  convenient  method 
for  computing  the  probabilities  of  conditions  for  different  types  of 
sort ies . 

This  report  presents  a  method  for  generating  these  probabilities. 

The  weapons  model  for  which  this  method  was  Intended  to  be  used  is  a 
linear  programming  model  written  by  Robert  Spelr  (Reference  1);  however, 
it  probably  can  be  used  for  many  other  weapons  allocation  models.  The 
method  consists  of  a  program  in  FORTRAN  IV  for  the  CDC  6600  which  converts 
data  given  by  the  Air  Weather  Service  (AWS)  for  a  given  location  into  data 
that  fit  the  launch  envelopes  of  various  weapon  systems.  The  location 
used  in  this  report  was  Berlin,  Germany. 


1.  Speir,  Robert  A.,  Another  Linear  Programming  Approach  to  USAF 
Conventional  Weapons  Planning.  AFATL-TR-73-201 ,  Eglln  AFB,  Florida,  1973. 
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SECTION  IT 


DESCRIPTION  OF  PROBLEM 


In  order  to  be  able  to  tell  a  weapons  allocation  model  whether  the 
weat'ier  is  good  enough  for  a  given  weapon  system  to  fly  its  sorties, 
there  must  be  a  reasonably  good  method  of  predicting  the  weather  at  the 
scene  of  the  battle.  Unfortunately,  methods  for  predicting  weather 
states  months  or  years  in  advance  do  not  exist;  hence,  only  probabilities 
can  be  dealt  with.  A  set  of  data  describing  the  probability  of  given 
weather  states  (in  terms  of  ceiling  and  visibility  conditions)  would  thus 
be  desirable.  Such  a  data  base  is  given  by  the  Air  Weather  Service.  For 
a  given  ceiling  level  C  and  visibility  level  V,  the  AWS  data  advance  the 
probability  that  in  a  given  3-hour  period  in  a  given  month  the  ceiling  will 
be  at  or  higher  than  C  and  the  visibility  will  be  at  or  higher  than  V. 
However,  this  may  not  be  compatible  with  the  requirements  of  the  allocation 
model.  The  model  may  require  ceiling  levels  or  visibility  levels  that  are 
not  given  in  the  AWS  tables,  and  the  conditions  for  using  certain  weapon 
systems  may  not  be  of  the  form  "celling  is  greater  Chan  or  equal  to  C  and 
visibility  is  greater  than  or  equal  to  V."  Hence,  a  method  for  converting 
Che  AWS  data  into  data  giving  the  probabilities  that  conditions  for  opera¬ 
tion  of  weapon  systems  are  favorable  would  be  desirable.  Such  a  method 
will  be  described.  It  consists  of  reading  Che  AWS  datA,  Interpolating  to 
get  the  desired  celling  and  visibility  levels,  converting  the  interpolated 
data  into  a  joint  probability  distribution,  averaging  the  data  over 
several  monthly  or  hourly  periods,  and  assembling  the  resulting  data  to 
give  probabilities  of  the  weather  states  given  in  the  linear  model. 
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Sl'.CTION  III 


MKTHOI) 


On  the  first  step,  tht  AV/S  liata  are  read  in.  These  data  consist  of 
a  set  C  =  of  celling  levels,  a  set  V  =  {Vl,**",Vn}  of  visl- 

hillty  levels,  and  a  cumulative  probability  distribution  F:CxV  [0,l] 
such  that  for  all  v  in  V  and  c  in  C,  the  number  F(c,v)  is  the  probability 
that  the  celling  is  -f  and  the  visibility  is  -V.  The  AWS  data  are  also 
dependent  on  the  month  m  and  the  hour  h,  so  that  F  is  really  a  function 
of  four  variables  c,  v,  m,  and  h. 

The  desired  input  to  the  linear  model  is  given  in  terms  of  condi¬ 
tions  on  the  ceiling  and  visibility  values,  such  as  "ceiling  between 
3,000  and  10,000  feet  and  visibility  between  3  and  5  miles."  Let  C'  be 
the  set  of  all  celling  levels  and  let  V'  be  the  set  of  all  visibility 
levels  referred  to  in  the  conditions.  In  general ,  C'  will  contain  some 
ceiling  levels  not  contained  in  C  and  V'  will  contain  some  visibility 
levels  not  contained  in  V.  The  data  must  be  expressed  in  terms  of  C' 
and  V,  and  to  do  this,  the  computer  program  Interpolates  between  the 
nearest  two  values  in  C  or  V  on  either  side  of  a  given  value  c'  in  C' 
or  v'  in  V'.  Two  interpolations  are  needed  -  one  for  ceilings  and  one 
for  visibilities. 

The  method  reported  here  uses  a  linear  interpolation  as  follows: 

Given  c'  in  C'  and  v'  in  V’,  suppose  that  c'  is  not  in  C  and  v*  is  not 
in  V.  The  level  c'  lies  between  the  two  nearest  levels  c~,  c"^  in  C, 
i.e.,  c"<c'<c+.  Similarly,  there  ai.e  the  two  nearest  levels  v+  and  v“ 
to  v'  in  V,  and  v"<v'<v‘*'.  (Assumption  is  made  that  zero  levels  are 
Included  in  C,  C' ,  V,  V',  and  that  the  highest  levels  in  C  and  V  exceed 
the  highest  levels  in  C'  and  V).  Define 

G'  (c,v')  =■  (F(c,v+)  -  F(c,v-))  +  F(c,v-)  (1) 

(G(c+,v)  -  G(c-,v))  +  G(c“,v)  (2) 

where  G;CxV'  -*•  [0,l]  is  an  intermediate  step  in  the  computation. 

Then  F'  is  considered  to  be  a  good  approximation  to  the  ceiling  and 
visibility  conditions  at  the  levels  contained  in  C'  and  V'.  This  con¬ 
stitutes  a  linear  Interpolation.  A  quadratic  approximation  would  not 
be  better,  as  demonstrated  by  the  following  example. 

In  one  actual  case  involving  Berlin  (visibility  6  miles,  October, 
0000-0200  hours)  the  following  probabilities  are  given. 


F’  (r',v') 
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IV  c 


Ceiling  I^evel  c 
(In  Feet) 


Probability  That  Visibility  ^6 
and  That  Ceiling  -c 


14,000 

12,000 

10,000 


0.266 

0.268 

0.284 


These  three  points  in  ceiling-probability  space  might  be  considered  in 
rder  to  obtain  the  probability  that  visibility  ^6  miles  and  that  ceiling 
13,000  feet.  Tlie  resulting  parabola  Is  given  by  the  equation 

P(l)  =  1.75x10"^  -  4.65xl0“‘’c  +  0.574 

It  is  required  by  the  laws  of  probability  that  the  requested  probability 
at  6  miles,  13,000  feet,  be  greater  than  the  probability  at  6  miles, 
14,000  feet.  However,  p(13,000)  =  0.2652  <  p(14,000)  =  0.266.  The  prob¬ 
lem  is  that  the  parabola  dips,  as  Figure  1  demonstrates.  This  sort  of 
situation  was  found  lo  occur  frequently  in  the  AWS  data.  Hence  a  linear, 
rather  than  a  quadratic,  approximation  was  used. 

Tf  v'  is  in  V,  then  the  value  G(c,v')  is  set  equal  to  F(c,v'). 
Similarly,  if  c'  is  in  C,  then  F'(c',v')  is  set  equal  to  C(c*,v*).  These 
two  steps  have  to  be  considered  separately,  since  applying  Equation  (1) 
or  (2)  will  result  in  division  by  zero. 


After  the  data  is  interpolated,  it  is  converted  into  a  joint  proba¬ 
bility  distribution.  To  do  this,  arrange  the  elements  of  C'  and  V'  in 
descending  order.  Then  v^  i  v2  -•••-  Vn  and  cl  c,„  (since  there  is 

no  longer  any  need  for  C  and  V,  the  primes  shall  be  dropped  on  the  v's 
and  c's).  This  implies  for  all  k  and  for  all  i  and  j  such  that  i  <  j, 
the  inequalities 

F  (ci,  vk)  <  F  (cj,  vk) 
and  F  (ck,  v^)  <  F  (ck,  Vj) 

hold.  Define 


and 


K  (Cj,  v^) 


F 


1,  Vj),  for 


i  =  m-1,***,  1, 

J  =  I,***,  n. 


(c^,  Vj) 


K  (c 


1’ 


K  (c^,  Vj-1).  For  j 
1 


n-1,*** ,  1 
!,•••,  m. 


where  K  is  an  intermediate  computational  function.  In  the  computer  program, 
F' ,  .1,  K  are  all  stored  in  the  same  array,  which  makes  It  necessary  to 
run  through  the  indices  backwards,  as  indicated  above.  For  the  remainder 
of  this  report,  F  will  denote  the  function  J  as  defined  above. 
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Figure  1.  Collocating  a  Parabola  Through  Three  Points  in 
Ceiling-Probability  Space 
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Tlie  next  step  is  to  average  the  data  over  a  set  of  hourly  periods  or 
over  a  set  of  months.  Two  alternatives  that  were  tried  were: 

(1)  Averaging  over  best  3  months  of  year  (or  over  the  worst). 

Best  was  defined  as  the  .)  months  with  the  highest  median  ceiling  (or 
highest  probability  of  no  ceiling,  in  case  the  median  celling  was  infinity). 
A  similar  definition  was  used  for  worst. 

(2)  Averaging  over  best  and  worst  3-hour  periods,  separately 
for  day  and  night.  A  3-hour  period  was  regarded  as  a  day  period  if  the 
sun  was  above  the  horizon  more  than  30  percent  of  the  time  during  the 
period;  otherwise,  the  period  was  a  night  period.  To  accomplish  this 
averaging,  separate  arrays  of  celling  and  visibility  probability  data 
for  each  of  the  month-hour  combinations  to  be  considered  were  read  in. 
Interpolated  and  converted  into  joint  distributions.  They  were  then 
averaged  together.  In  the  results  given  in  this  report,  no  averaging 
took  place. 

The  final  step  consisted  of  totaling  the  present  set  of  data  so 
that  probabilities  for  the  conditions  required  by  the  linear  model  could 
be  calculated.  These  conditions  can  be  expressed  by  giving  a  partition 
of  the  set  C'xV'  : 

C  X  V  =  L,  U  L„  u***u  L  , 

12  r 

where  the  unions  are  disjoint.  The  L's  can  be  thought  of  as  weather 
states.  For  example,  if  C'  *=  {lOOOO,  6000,  AOOO,  0}  and  V'  ■  {10,  6, 

4,  0} ,  then  a  possible  partition  is  as  follows: 

*  {  (10000,  10),  (10000,  6),  (10000,  4),  (10000,  0)  } 

=  {  (6000,  6),  (6000,  4),  (4000,  6),  (4000,  4)  } 

=  {  (6000,  10),  (4000,  10)  } 

-  {  (6000,  0),  (4000,  0),  (0,  10),  (0,  6),  (0,  4),  (0,  0)  } 

The  interpretation  of  and  is  as  follows: 

Lj^  Ceiling  -10,000  feet 

L2  Ceiling  between  4,000  and  10,000  feet  and  visibility 
between  4  and  10  miles 

L-  Ceiling  between  4,000  and  10,000  feet  and  visibility 
-10  miles 

L,  Ceiling  14,000  feet  or  visibility  l4  miles  and  celling 
^  110,000  feet 
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As  the  above  descriptions  show,  the  Lj  can  be  described  by  specifying 
for  each  ceiling  level  the  vlsiliillty  levels  contained  in  it.  The 
algoritlim  then  adds  up  the  probabilities  in  each  of  the  hj^'s  to  obtain 
the  probability  that  each  weather  state  will  hold,  in  symbols, 

p  (L  )  =  1  F(c,v)  for  i  =!,•••,  r 

(c,v)cl_ 

These  probabilities  can  then  be  used  as  data  for  the  linear  programming 
model  of  Spelr  (Reference  1). 
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SECTION  IV 


INPUT,  OUTPUT  ANT)  EXAMPLE 

To  demonstrate  the  execution  of  this  program,  consider  the  following 
problem:  Suppose  a  certain  military  operation  is  to  be  undertaken  between 

OOnO  and  0300  on  5  October  1980  near  Berlin,  Germany.  To  determine  how 
the  operation  is  to  take  place,  it  is  desired  to  find  the  probability  of 
the  following  weather  states: 

(1)  Ceiling  -  344  feet  and  Visibility  -  2.56  miles 

(2)  344  feet  1  Ceiling  ^  1784  feet  and  Visibility  ^  2.56  miles 

(3)  1784  feet  1  Celling  -  3033  feet  and  Visibility  i  2.56  miles 

(4)  Celling  ^  3033  feet  and  2.56  miles  £  Visibility  £  5.12  miles 

(5)  3033  feet  £  Celling  £  9874  feet  and  Visibility  ^  5.12  miles 

(6)  Celling  ~  9874  feet  and  Visibility  -  5.12  miles. 

These  are  mutually  exclusive.  The  data  to  be  used  in  determining  the 
above  probabilities  is  the  AWS  data  for  October,  0000-0300  hours. 

The  input  to  be  provided  to  the  computer  is  the  following: 

1.  Card  1  (12,  II,  412) 

MTV  Month 

IHR  Hour  Period  (0000-0300  is  1,  0300-0600  is  2,  etc.) 

NC  Total  number  of  ceiling  levels  referred  to  in 

description  of  weather  states 
NV  Total  number  of  visibility  levels  referred  to  in 

description  of  weather  states 
IT  Total  number  of  weather  state  descriptor  cards 

(see  below) 

NWC  Number  of  weather  states. 

2.  Card  2  (13F6.1) 

The  NC  ceiling  levels  (in  feet)  in  order  from  highest  to 
lowest  (XC(1),  1*1, N). 

3.  Card  3  (13F6.1) 

The  NV  visibility  levels  (in  miles)  in  order  from  highest 
to  lowest  (XV(I),  1*1, N). 
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A.  Ciirds  A  to  IT  +  3  (AlAj 

TT  cards  describing  the  weather  states  as  follows; 

NN(I)  Number  of  weather  state 

NCN(l)  CelHng  zone  (XC(I)llevel^(T-l))  for  which  this 
card  gives  the  visibility  zones  corresponding  to 
the  weather  state  I. 

NVA(I)  (The  visibility  level  considered  in  this 
NVR(I)  (card  is  NVB(I)llevel£NVA(I) . 

As  an  example  for  this  report  the  weather  states  can  be  Indicated  by  the 
following: 


Visibility  Levels  (Miles) 


Celling  Levels 
(Feet) 


12  3 

5.12+  2.56-5.12  0-2.56 


987A+ 

59A2  -  987A 
3033  -  59A2 
178A  -  3033 
3AA  -  178A 
0  -  3AA 


6 

5 

5 

3 

2 

1 


A 

A 

A 

3 

2 

1 


1 

1 

1 

1 

1 

1 


Thus,  to  describe  weather  state  5,  the  following  cards  are  used: 

5  2  11 
5  3  11 

That  is,  weather  state  5  consists  of  ceiling  level  2(59A2-987A  feet)  at 
visibility  level  1  (5.12+  miles)  and  ceiling  level  3(3033-59A2  feet)  at 
visibility  level  1  (5.12+  miles).  Thus,  the  complete  set  of  data  for 
this  problem  is  as  follows: 

101_6_31A_6 

987A.  59A2.  3033.  178A.  3AA. 


5.12 

2 

.56 

_ ^0. 

1 

1 

3 

1 

"2 

3 

3 

1 

■3 

3 

3 

1 

"a 

3 

3 

1 

"5 

3 

3 

1 

"6 

1 

3 

2 

”5 

1 

2 

3 

”a 

1 

2 

A 

”1 

2 

2 

0. 


A _ 2  _2 _ 2 

A _ 3'~  2 _ 2 

5 _ 2 _ 1 _ 1 

5 _ 3  _ 1 

'6  l"  1  1 


The  original  weather  data  should  be  on  a  TAPE  A  in  binary  format,  giving 
first  the  month  number,  then  the  hour  period  number,  then  a  series  of 
data  representing  a  32x16  array  containing  the  AWS  data  for  the  given 
month  and  hour  period.  The  tape  should  give  AWS  data  for  all  month-hour 
period  combinations.  The  program  advances  TAPE  A  until  the  correct 
month-hour  period  is  found  according  to  the  first  data  card.  The  program 
was  run  on  a  CDC  6600  with  Scope  3.3;  on  another  computer  system  some 
modifications  of  the  program  may  have  to  be  made. 

The  output  has  four  parts: 

1.  Joint  probability  distribution  of  ceiling  versus 
visibility  for  the  ceiling  and  visibility  data 
given. 

2.  Conditional  probaUllty  of  the  ceiling  levels  given 
a  certain  visibility  level. 

3.  Conditional  probao Mty  of  the  visibility  levels 
given  a  cert  tin  ceiling  level. 

A.  Probabilitiw.-  of  the  given  weather  states. 

Note  that  in  the  output  the  figure  5942  as  a  celling  level  means 
celling  between  59A2  feet  and  the  next  highest  celling  level,  i.e., 

9874  feet.  A  listing  of  the  FORTRAN  IV  program  followed  by  results 
for  this  example  are  included  in  Appendix  A. 
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SKCTTON  V 


CONCLUSIONS 


The  purpose  of  this  report  was  to  provide  a  method  by  which  weather 
data  provided  by  AWS  ran  be  converted  Into  a  form  usable  in  weapons 
allocations  programs,  with  Speir's  program  (Reference  1)  specifically 
In  mind.  The  methodology  in  the  report  is  more  important  than  the 
specific  details.  Although  the  original  AWS  data  tapes  have  been 
purged  and  the  data  available  only  in  table  form,  the  program  or  a 
suitable  modification  of  it  will  provide  weather  state  data  most  of 
the  time  when  it  is  needed. 


11 

(The  reverse  of  this  page  I*;  blank) 


APPENDIX  A 


FORTRAN  IV  WEATHER  STATE  PROGRAM 
AND 

RESULTS  OF  WEATHER  STATE 
EXAMPLE 
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I 


INITIAL  DISTRIBUTION 


HQ  USAF/SAMI  1 
AFIS/INTA  1 
ASD/ENFEA  1 
00-ALC/MMWM  2 
TAWC/TRADOCLO  1 
AFATL/DL  1 
AFATL/DLODL  9 
AFATL/DLODR  1 
AFATL/DLY  I 
AFATL/DLYD  1 
AFATL/DLYW  1 
TRADOC  (ADTC-DO)  L 
DDC  2 
AUL  (AUL-LSE-70-239)  1 
TAC/INA  1 
TAC/DRA  1 
ASD/XRP  1 
PACAF/DOO  1 
US  Army  TRADOC  1 
USAFE/DOQ  1 
ASD/ENYS  1 
AMXLE-PM  I 
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